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ABSTRACT 

Obtaining the magnetic properties of non-resolved structures in the solar photosphere is always 
challenging and problems arise because the inversion is carried out through the numerical minimization 
of a merit function that depends on the proposed model. We investigate the reliability of inversions in 
which the stray-light contamination is obtained from the same observations as a local average. In this 
case, we show that it is fundamental to include the covariance between the observed Stokes profiles and 
the stray-light contamination. The ensuing modified merit function of the inversion process penalizes 
large stray-light contaminations simply because of the presence of positive correlations between the 
observables and the stray-light, fundamentally produced by spatially variable systematics. We caution 
that using the wrong merit function, artificially large stray-light contaminations might be inferred. 
Since this effect disappears if the stray-light contamination is obtained as an average over the full 
field-of-view, we recommend to take into account stray-light contamination using a global approach. 
Subject headings: methods: data analysis, statistical — techniques: polarimetric — Sun: photosphere 



1. INTRODUCTION 

The quantitative investigation of the magnetism of 
structures in the solar atmosphere is done through the 
analysis of the observed Stokes profiles. Particularly dif- 
ficult is to infer the properties of weakly magnetized 
regions of the solar surface such as the quiet internet- 
work. The reason is that the observed polarization sig- 
nals often stay at the detection limit of modern spectro- 
polarimeters because the spatial resolution is still not 
high enough to resolve the smallest magnetic structures. 
Therefore, the magnetic properties are inferred assuming 
unresolved structures, which is model dependent. 

Due to lack of information, it is customary to apply 
relatively simple models to explain the observations. The 
model parameters are estimated using a least-squares 
(or, equivalently, maximum- likelihood) approach (e.g., 
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many others) or using a fully Bayesian app roach 
(jAsensio Ramos et all l2007t lAsensio Ramos! [2009) . For 
the inversion of unresolved structures, the idea of using 
several components that contribute to the observed 
Sto kes profiles has bee n arou nd since its introduction 
by iLites k, Skumanichl (|1990l) . In its simplest form, 
the method consists in using a filling factor a of the 
pixel that accounts for a non-magnetized component 
(cither stray-light or a pure non-magnetic plasma) and 
the remaining 1 — a fraction of the pixel is filled by a 
magnetic component. 

This paper points out that, in case the stray-light con- 
tamination is obtained directly from the observations as 
an average over a given field-of-view, some modifications 
are necessary in the inversion procedure to take into ac- 
count the eventual presence of correlation between the 
observed Stokes profiles and the stray-light. This modi- 
fication has the remarkable property of penalizing large 
stray-light contaminations. 



2. DISCUSSION 

In order to obtain information on the magnetic and 
thermodynamical properties of the solar atmosphere, one 
proposes a model atmosphere that depends on a set of 
parameters 6 (e.g., the temperature T at one or sev- 
eral heights, hydrogen density, magnetic field strength 
and inclination, etc). This model is used to synthesize 
the Stokes vector O mod (a;, y, A; 6) for an arbitrary num- 
ber of spectral lines at given spatial position (x,y) and 
wavelength A. Expressing the observed Stokes profiles as 
O obs (a;, y, A), it is customary to obtain the "best" param- 
eters as those m inimizing the following merit function 
(e.g. JPress et al.lll986( l: 
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where the sum over j is extended to all N observed wave- 
length points. In this equation, of (a;, y, Xj) represents 
the variance of the numerator, for each wavelength Xj, 
Stokes parameter i = I, Q, U, and V, and at each lo- 
cation (x,y), due to possible uncertainties in the obser- 
vations (measurement errors and noise). This is differ- 
ent from "real" variations due to space-time fluctuations 
(e.g., intensity contrast due to granule-intergranule fluc- 
tuations). Usually, the model is known with certainty 
— as when we fit a Gaussian profile to an observed spec- 
tral line — , and of (x, y, Xj) is just the noise variance a\ 
of the data. This noise variance can be estimated from 
the observations, ideally, by taking several observations 
of the same object under identical observational condi- 
tions. Unfortunately, this is often just not possible. It is 
then customary in spectroscopic observations to select a 
continuum window and, assuming that it should be spec- 
trally flat, all fluctuations are due to noise. The flatness 
assumption is usually not fulfilled due to systematic ef- 
fects and the estimated variance might be larger than the 
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Figure 1. Value of the variance in units of as defined in Eq. 
for different values of M and for four different values of p, the 
covariance between the observed profile and the stray-light. Note 
that p depends on M but, for simplicity, we have assumed that 
Cov [0° hB (x, y), 0° hB (x' , j/')] is constant, so that p is independent 
of M. Although not shown to avoid crowding, when the covariance 
is maximum (p = 1), values of the stray-light contamination above 
~l/2 are strictly forbidden. The exact value of this limit can be 
obtained from Eq. J8T . 

one associated with random noise effectsQ 

However, it is also possible that the model contribute 
to the variance uncertainty by incorporating explicitly 
some observable — for example, if we fit a Gaussian pro- 
file to an observed spectral line, but keeping the equiva- 
lent width of the observed line profile. We study here an 
example of this type that arises naturally when consid- 
ering stray light in spectropolarimetric observations. 

The proper interpretation and inversion of spectropo- 
larimetric observations with spatial resolution requires 
some treatment of stray light. "Stray light" refers to the 
unavoidable spread of light from different regions on the 
source. It is an effect produced by the extended tails of 
the telescope's point spread function, which causes that 
a significant fraction of the photons detected in one pixel 
results from regions outside the Airy disk. As a conse- 
quence, there is always some amount of contamination 
from regions other than the region that was supposed to 
be imaged by the instrument. This problem has become 
especially pressing recently, with the advent of high- 
resolution space-b orne spectropola rimeters like the SP 
(jLites et al.1 120011 ) aboard Hinode (|Kosugi et all 120071 ). 
To deal with this problem when carrying out inversions 
of Stokes profiles from quiet regions of the solar pho- 
tosphere, the following model has been proposed (e . 



i psed (e -g., 
1 l27)TO f) : 



ILites fc Skumanichl[T990l : lO rozco Suarez et al 

O mod (x, y; 9) = aT>(x, y) + (1 - a)S(x, y; 6), (2) 

where S(x,y,6) are the Stokes profiles emerging from a 
model atmosphere representative for a magnetized region 
occupying a fraction 1 — a of the pixel (from now on, 
we drop the dependence of the variables on A). The 
term T)(x,y) is a stray-light contamination profile that 
is obtained from a local average of the M pixels around 
the pixel of interest (of the order of 1 arcsec 2 around the 

1 We shall not consider here the expected variation of er^ with 
wavelength across strong spectral lines. This is a very interesting 
problem that will be considered elsewhere. 



pixel of interest for Hinode observations) : 

T>(x,y)=M- 1 J2 °° b V>2/), (3) 

(x\y')en 

where fl is the set of M pixels that is considered to be 
affecting the pixel of interest. The parameter a quantifies 
the amount of stray-light contamination. In this case, the 
denominator of the x 2 function in Eq. (TT]) is: 

oj(x, y) = Var [aD z (x, y) + (1 - a)Si{x, y) - 0° hs (x, y) 

(4) 

The presence of Di(x,y), which depends on the ob- 
servations 0° hs (x\ y'), produces that the variance is not 
just the noise variance because a certain degree of corre- 
lation might exist between 0° bs (x, y) and Di(x,y). The 
reason is that both quantities are affected by noise and 
systematic effects and they are the result of a complex re- 
duction process that might introduce some correlation. 
A straightforward calculation allows us to simplify the 
variance to: 

a 2 (x,y) = (l + 50 a 2 n - 2aCov [0° bs (x, y), D t {x, y)} . 

It is of interest to define p as the covariance measured 
in units of the noise variance: 



pa 2 n = Cov[0° hs (x,y),D l (x,y)] 

' J2 Cov[0? s ( X ,y),0?*(x',y')], (6) 



M 
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which explicitly shows the dependence of Di(x, y) on the 
observed Stokes profiles and also indicates that p behaves 
like an average correlation coefficient. Plugging this ex- 
pression into Eq. ([S]) results in: 



o-t(x,y)= [l + — -2ap 



(7) 



This expression for the variance has several interest- 
ing properties. It depends on the stray-light contami- 
nation coefficient a. The term in a 2 accounts for the 
inclusion of additional noise variance coming from the 
stray-light profile. Its influence is heavily reduced when 
averaging over many pixels, decreasing in proportion to 
the number of pixels added. The linear term depends on 
the covariance between the observed Stokes profile at a 
given wavelength and the average profile that is consid- 
ered as stray-light contamination. This term also goes to 
zero when many pixels are considered for the averaging, 
so we recommend to obtain the stray-light contamina- 
tion profile as an average over the full field-of-view. A 
non-zero contribution can be produced by any spatially 
variable systematic effect present in the final focal plane 
produced by the camera (flatfield systematics or fringes) 
or by any optics before (fringes). Additionally, correc- 
tions carried out during data reduction can introduce 
correlations between surrounding pixels. Finally, data 
compression like t he JPEG compression used by Hinode 
(|Lites et al.l 120021 ) to optimize telemetry can also intro- 
duce correlation artifacts. As a consequence of this de- 
pendence on a, the \ 2 merit function to be optimized 
to get the maximum-likelihood parameters 6 is different 
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Figure 2. Spatial variation of the variance (upper left panel) and covariance between the local profile and the stray-l ight contam i natio n 
using a local environment of 1.5" around each pixel (upper right panel). The Hinode observations are those used bv lLites et a!. , (2008). 
This figure is representative of a window in the continuum and it is expected to change if other points in the spectral direction are used 
because the systematic effects might change. The lower panels show the histogram of each image. Note that \p\ > 1 are found because 
of the small sample with which it is calculated, thus leading to statistical fluctuations in which the covariance and the variance are not 
strictly compatible. 

from the one commonly used in the literatur^E We shall 
see shortly that this can affect the results significantly. 

Although —1 < p < 1, strong anticorrelations (p < 0) 
can be effectively discarded in our case. Figure [T] shows 
how of/ (7^ [see Eq. ((7|)] varies with M and p, assuming 
that Cov [0° hs (x, y), 0° hs (x', y')] is constant on Cl. Tak- 
ing M = 24 is roughly equivalent to consider a neigh- 
borhood ~1.5" around the pixel of interest at Hinode's 
spatial sampling. As seen in Fig. [1] the dependence of 
the variance on M is relatively weak and the parabolic 
shape of the curves tends to a straight line as M in- 
creases. If p > 0, stray-light contaminations above the 
following limit arc forbidden: 



= pM - ^{pMf 



(8) 



which has been obtained as the solution of af(x,y) = 
0. This quantity rapidly tends to (2p) _1 for increasing 
values of M and is smaller than 1, whenever p > 1/2 
and M > (2p - If p < 0, the behavior is the 

opposite and large values of the stray-light contamination 

2 Note that this modifies the likelihood function used by 
Ascnsio Ramos (2009) in the Bayesian framework too. 



are favored. 

It is important to estimate the covariance between 
stray-light contamination and observed profi les. To this 
end, w e use the quiet Sun map observed by iLites et al.l 
( 2008) to compute the stray-light contamination in a typ- 
ical Hinode observation as an average over a window of 
~1.5" around every pixel of interest. The upper left 
panel of Fig. [5] shows the variance of the observed Stokes 
/ continuum normalized to the average variance of the 
whole map. This is equivalent to the noise variance in 
case the continuum is assumed to be spectrally flat. Our 
experience is that the average value of the standard de- 
viation for Stokes / is a factor 3-4 larger tha n the noise 
estima tion in Stokes Q, U and V made by Lit es et all 
(2008). We assign this difference to the unavoidable pres- 
ence of systematic effects in Stokes / due to flatfielding 
procedures. Since they change from pixel to pixel, they 
can be absorbed as part of the noise, though probably 
not normally distributed. Note the presence of conspic- 
uous horizontal stripes that show the presence of pixels 
along the slit and at the points of the camera associated 
with continuum wavelengths which present a somewhat 
higher variance (higher systematic effects). The upper 
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right panel of Fig. [3] shows the value of p for the whole 
map, which has been calculated using a small window in 
the continuum from 6302.915 A to 6303.002 A. Strictly 
speaking, the covariance should have been calculated us- 
ing many realizations of the measurement process with 
the underlying solar profile fixed. Assuming some kind 
of ergodicity, we estimated the covariance using a wave- 
length window in the continuum. For consistency, we 
also estimated the covariance al ong the time using the 
time series of (jLites et al.l [2008) , with very similar re- 
sults. Again, the presence of horizontal stripes is clear, 
pointing to the systematic character of such defects. On 
average, p ~ 0.45, although the distribution of covari- 
ances is clearly heavy-tailed, with the presence of large 
values of p much more frequently than in the case of a 
Gaussian distribution. The lower panels of Fig. [5]present 
histograms of the maps. A characteristic is that it is 
possible to find unphysical \p\ > 1, a consequence of the 
small sample with which this quantity is computed. 

It is customary in standard inversion codes to intro- 
duce different weights for each Stokes parameter so that 
Stokes / does not dominate the merit function. One 
may think that this alleviates the effect of covariance 
on the inferred parameters. However, the inclusion of 
weights does not mimic properly the form of x 2 because 
the variance term in Eq. ((TJ) depends explicitly on one 
of the model parameters, although the intuition is par- 
tially right. Therefore, one would expect differences with 
respect to the parameters obtained with a standard x 2 - 
This is indeed the case, as shown in Tab. [TJ This table 
shows the results of a least-squares fit to the Stokes pro- 
files shown in Fig. [3] for different values of p. The model 
proposed is that of Eq. where the magnetic com- 
ponent is a Milne-Eddington atmosphere and the stray- 
light contamination is obtained from a local average of 
M = 24 pixels around the one of interest. The table 
presents the magnetic field strength B and inclination 
8b, together with the stray-light contamination a that 
minimizes the x 2 for different values of p and the value of 
the reduced \ 2 a t the minimum. A striking effect of the 
inclusion of stray-light is that the value of a decreases as 
soon as p increases, a direct consequence of the presence 
of positive correlations between the observed profile and 
the stray-light. Therefore, the magnetic field strength 
has to be adjusted accordingly in order to maintain the 
magnetic flux de nsity, much in the direction of wh at has 
been explored bv lMartmez Gonzalez et all (|2006l) . Note 
that the cases with p > 0.5 are not so common in the 
observations. As a consequence, since such large values 
are not representative of the sample profile, the fits are 
not good and they should be taken as indicative of what 
would happen in such a limiting case. Finally, it is im- 
portant to point out that this behavior would be reduced 
if the field is strong enough to produce a substantial Zee- 
man splitting because the magnetic field strength can be 
inferred directly from the splitting. 

3. CONCLUSIONS 

We have demonstrated that the inversion of Stokes pro- 
files with models that include stray-light contaminations 
obtained from the same observations has to be carried 
out with care. In general, the x 2 merit function is mod- 
ified and has to include the effect of the variance of the 
stray-light profile and the covariance between such pro- 
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Figure 3. Observed Stokes profiles (black solid line) of a pixel in 
the field-of- view observed by ILites et al.l j2008j) . Inversions with 
p = are shown in red, p = 0.5 in blue, p = 0.7 in green and p = 1 
in magenta. All of them give reasonably good results but the model 
parameters change dramatically as shown in Tab. [T] Note that the 
cases of large p are not realistic in this situation since we adopt a 
covariance that is too large with respect to the correct one. 



Table 1 

Maximum-likelihood parameters 
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0.91 
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0.88 
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0.7 


114 


176 


0.46 


11.3 


1.0 


62 


169 


0.02 


13.8 



file and the one being inverted. The most important 
consequence is that the noise variance now depends on 
the stray-light contamination parameter, a. 

A first contribution adds quadratically with the noise 
variance and is inversely proportional to the number of 
pixels that contribute to the stray-light profile. When 
the number of pixels is large enough, this term turns out 
to be almost negligible. The fundamental reason is that 
the signal-to- noise ratio of the stray-light profile increases 
and becomes much larger than that of the local profile, 
which dominates then the variance. If the stray-light 
contamination profile is obtained as an average over the 
full ficld-of-view, the effect of correlation diminishes con- 
siderably and the modified x 2 converges to the standard 
X 2 used in the past. 

A second contribution adds linearly and is much more 
delicate. It accounts for all the systematic effects 
(fringes, flatfield systematics, data reduction effects, etc.) 
that plague the observations. The presence of systemat- 
ics that are well above the noise level make both profiles 
vary similarly. The covariance term takes this into ac- 
count by disfavoring models with large stray-light con- 
taminations. In the absence of this term, any inversion 
method prefers to use a large stray-light contamination 
because the profiles become much more similar using the 
standard and incorrect \ 2 ■ A simulated exercise have 
shown that neglecting the presence of this covariance 
might lead to artificially strong magnetic field strengths 
and large stray-light contaminations. 
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